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Light propagation in a photonic crystal infiltrated with polarizable molecules is considered. 
We demonstrate that the interplay between the spatial dispersion caused by Bragg diffraction and 
polaritonic frequency dispersion gives rise to novel propagating excitations, or braggoritons, with 
intragap frequencies. We derive the braggoriton dispersion relation and show that it is governed by 
two parameters, namely, the strength of light-matter interaction and detuning between the Bragg 
frequency and that of the infiltrated molecules. We also study defect-induced states when the 
photonic band gap is divided into two subgaps by the braggoritonic branches and find that each 
defect creates two intragap localized states inside each subgap. 



I. Introduction 

Photonic crystals and- in particular photonic band 
gap (PEGU. materialsEra have recently attracted much 
attentionEm due to their rich physics and possible appli- 
cations. In these systems the dielectric function is period- 
ically modulated and, as a result, their optical properties 
are dominated by light diffraction effects. When Bragg 
diffraction conditions are met then light scattering is very 
strong, so that within certain frequency intervals near the 
resonances light propagation is inhibited. 

Sinc£pthe subject of photonic crystals was intro- 
ducedJU'El one of the main goals of photonic band- 
structure calculations was to engineer structures with a 
complete band gap, i.e. with no propagating solutions 
of Maxwell's equations within a certain forbidden gap. 
The pursuit of this goal has generated a stream of stud- 
ies that are too numerous to be cited here; early works 
are reviewed in Refs. ^|, ||. Here we only mention that a 
complete bandrgap in two dimensions (2D) was theoreti- 
cally predictedU'B and experimentally demonstrated^ for 
an array of dielectric rods. In the quest for a structure 
having a complete PBG in- three dimensions (3D), the 
diamond lattice was showro to be morepjpromising than 
a simple face centered cubic (fee) latticeEm. 

The frequency gap in the photonic spectrum sets a 
stage for a number of physical effects. The prime ef- 
fect, namely the inhibition of spontaneous emission for 
an emitter with transition frequency within .-.the, gap, 
was already suggested in the pioneering works .tlclll3 Fur- 
thermore, since light cannot leave the emitting atom, a 
coupled atom-field in-gap state is formed, in which the 
atomic level is "dressed" by its own exponentially local- 
ized radiation field. Blill It was also demonstrated that 
although a single photon cannot propagate inside the 
gap, nevertheless a non-linear medium embedded inside 
the photonic crystal gives-, rise to multi-photon bound 
statesLL3, or gap solitona!3 that result in self-induced 
transparency. Yet another consequence of PBG is the 
modification of cooperative emission with frequency close 
to the band edge. In particular PBG was shown to 



change the rate o f a up erradiant emission from an en- 
semble of emitters. EiJli-3 Lastly, PBG structures facilitate 
strong Anderson localization of photonsta because the 
sharp density of states within the gap spectral range n&* 
cessitates a reinterpretation of the loffe-Regel criterion.tl 

PBG structures with a defect constitute a separate 
area of study initiated by the classical works in Refs. 
||, |l6|. These structures are important since the defects 
cause localized intragap states. For these states, the PBG 
sample acts as a resonator with a very high quality fac- 
tor. This property was recently used for designing a low- 
threshold PBG defect-mode laser .O 

Another class of materials with a forbidden gap 
for light propagation are spatially homogeneous, but 
frequency-dispersive media. The energy gap in these sys- 
tems has a polaritonic origin, i.e. it is formed due to 
the interaction of light with the medium polarization.Ea 
This energy gap can be viewed as the result of anticross- 
ing between the photonic and excitonic dispersion rela- 
tion branches. Some non-trivial manifestations of the 
polaritonic gap were recently explored in Refs. |l^, In 
these papers a general model of two-level systems inter- 
acting with elementary electromagnetic excitations with 
a gap in the spectrum was solved by means of the Bcthc 
ansatz technique. Within t hi s m odel a very rich exci- 
tation spectrum was found,Ejn3 consisting of ordinary 
solitons, single-particle impurity bound states and mas- 
sive pairs of confined gap excitations and their bound 
complexes — dissipationless quantum gap solitons. 

Most of the available photonic crystals nowadays how- 
ever have incomplete PBGs; this means that light prop- 
agation is forbidden only along certain directions in- 
side the crystal. A prominent example are opalsp— rep- 
resenting self-assembled monodispersed silica ballsEll ar- 
ranged in a fee type lattice. Although opals have only 
incomplete PBG, the voids between the balls can be 
infiltrated by various media, which brings about non- 
trivial physics. In particular, the medium may con- 
tain polarizable molecules. Infiltrated opal with polariz- 
able molecules combines therefore polaritonic and Bragg- 
diffractive properties. Obviously, both effects coexist in- 
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dependently when the Bragg (ui = u>b) and polaritonic 
(lu = lo'y) resonances are well separated in frequency. 
A completely different situation occurs when lub ~ wt- 
This may be easily achieved in infiltrated opals that gives 
rise to a peculiar interplay between various frequency dis- 
persions. This interplay is the subject of the present pa- 
per. 

Our most important finding pertains to the case when 
the polaritonic gap of the polarized molecules infiltrat- 
ing the opal lies within the opal PBG. We demonstrate 
that such an overlap gives rise to novel massive propagat- 
ing excitations having frequencies inside the Bragg gap, 
that we dubbed here braggoritons. In other words, the 
Bragg gap splits into two sub-gaps, so that the braggori- 
tonic branches are isolated from the rest of the spectrum. 
We found that the braggoriton dispersion relation is very 
sensitive to the frequency detuning between wb and ut 
and to the relative width of the polaritonic gap (or, al- 
ternatively Rabi frequency) and the Bragg gap. 

The principal assumption we adopt here is that the 
Bragg gap, Alo-q, is narrow compared to wb; this is actu- 
ally the case in opals. The small value of Acjb/^b (<C 1) 
enables us then to obtain analytical results. In addition 
we also study the phase slip related intragap defect states 
for wb ~ wt- In the absence of polaritonic effect, the un- 
derlying physics of the defect-induced intragap states was 
already discussed in the original PBG paperEl. An anal- 
ogy was drawn between a defect state and a localized 
mode in a distributed feedback resonator, that originates 
from a phase slip. We extend this picture to incorporate 
polarizable medium and show that when the Bragg gap 
splits into two sub-gaps, then an existing phase slip gives 
rise to two localized states with frequencies within each 
of the sub-gaps. 

Our paper is organized as follows: In order to introduce 
the notations, we separately review in Sec. II the deriva- 
tion of the PBG and polaritonic spectra using the second 
quantization representation. In Sec. Ill we consider the 
combined Hamiltonian in the second quantization repre- 
sentation and diagonalize it by a unitary transformation. 
This yields the dispersion relations for the two excita- 
tions outside the gap, or Bloch-like waves, and the two 
intragap branches, or braggoritons. The properties of 
these novel excitations are analyzed in Sec. IV. We use 
them in Sec. V to determine the intragap frequencies of 
the defect-induced localized states. Concluding remarks 
are presented in Sec. VI. 

II. Second quantized PBG and po- 
laritonic Hamiltonians 

The Hamiltonian TL of the system under study is the sum 
of three terms 

H = H p h + W m + Hm-ph- (1) 



The first term, 7i p h, describes the photons in a pho- 
tonic crystal. The second term, Ti m , is the Hamiltonian 
of the polarizable medium; Ti. m - P h describes the photon- 
medium coupling. In this Section we review two limiting 
cases: (?) no polarizable medium (7i m = 0), and (u) no 
modulation of the dielectric constant. 

A. Incomplete PBG 

The general form of the Hamiltonian 7i p h is 

«ph = -tj dr [ £ ^)£ 2 + H 2 ] , (2) 

where E and H are respectively the electric and mag- 
netic fields. For a constant .dielectric function, e(r) = Eq, 
the second quantized formEa of the Hamiltonian (2) re- 
duces to a sum over oscillators representing plane waves 
with frequencies ujk = ck/^/eo, where k is the wave vec- 
tor. Modulation of e(r) causes light diffraction, so that 
the plane wave solutions are no longer the correct eigen- 
functions of the Hamiltonian (2). Below we consider 
a photonic crystal with an incomplete PBG along the 
z axis. This situation can be approximated by a one- 
dimensional modulation of e(r) along the z direction. If 
the Bragg gap is relatively narrow as assumed above, then 
it may be sufficient to consider only the first harmonics 
in e{z): 

e(z) — e + 6e cos(crz + 0). (3) 

Here Se (^C Eq) is the modulation amplitude, a = 2-n/d, 
where d is the modulation period and (f> is the dielec- 
tric modulation phase. We assume for simplicity that 
the electromagnetic field propagates along the z-direction 
and is homogeneous in the xy plane. In this case light po- 
larization is irrelevant. Generalization to arbitrary prop- 
agation direction is straightforward. The Fourier compo- 
nents of e(z) in (3) couple the original photon oscillators 
with momenta k and k ± a. These coupled oscillators 
form an infinite series that is constructed by successive 
addition (subtraction) of a. However, if Se <C £o and 
the wavevectors domain is restricted to the vicinity of 
the first Bragg resonance at k w cr/2, then the Hamil- 
tonian (2) can be truncated. In this case only the cou- 
pling to the near-resonance backscattered photons with 
momenta (a — k) cr/2 must be retained, so that the 
Hamiltonian (1) takes the form 

W p h = 

^2 {w(<?)a+ (q)a^(q) + w{-q)at(-q)a^{-q) 
q 

+ n B [e* a+(q)a^(-q) + e** a+ (-g)o^(g)] } . (4) 

Here, we introduced the notations: q = fc — a/2, a^(q) — 
a-k and a^(—q) = au- a for k w a/2, where dk is the 
usual photon annihilation operator. In the notations in- 
troduced in (4), the frequencies of the photonic branches 
are given by: 
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c(q + a/2) ( 2q 
uj(q) = — = wb H 



(5) 



where wb = ca /(2y/eo) is the Bragg frequency. We de- 
fine the coupling constant, Ob, as the half- width of the 
Bragg gap, i.e. Ob = ^Awb- It can be shown that Ob is 
related to the amplitude of the dielectric function mod- 
ulation: Ob = WB<5e/(2eo)- The summation in Eq. (4) is 
performed over the k domain \q\ -C cr/2. 

It is straightforward to diagonalize the Hamiltonian in 
Eq. (4) with the use of the following unitary transforma- 
tion 



(q) = cos 6 A (q) + sin 6 e** fa (q) , 
a_(-g) =-sinl? e _i * pi(q) + cos6> p 2 (q), 



where 



cos26> 



U)(q) - uj(-q) 



y/[w(q)-Lj(-q)] 2 + 4n 



(6) 



(7) 



The new operators B\ and B2 describe the creation (an- 
nihilation) of pairs of Bloch waves, that consist of forward 
and backscattered photons near the Bragg frequency. 
The diagonalized Hamiltonian (4) takes the form 

^ P h - E w b 1} (?) Pi (?) A (?) + 4 2) (?) , (8) 



where the dispersion relations of the two photonic 
branches are given by: 



,(1,2) 



(?) 



40? 



0|. 



(9) 



As mentioned above, the width of the PBG from Eq. (9) 
is 20b- 

B. Polarizable medium 

The Hamiltonians 7i m of polarizable medium and 7i m _ p h 
of light-polarization interaction in Eq. (1) can be written 
in the second quantization form as: 



k 



(10) 



assume here to be dispersionless having frequency wt; 
Op denotes the light-medium coupling strength that is 
proportional to the Rabi frequency. In the absence of 
the Bragg scattering term (Ob = 0) the complete Hamil- 
tonian (l)reduces to the conventional polaritonic Hamil- 
tonian 7i po i, given by: 



k 



k + at b k 



(11) 



with eigenstates representing the mixture of light and 
medium excitations 



Sfe = cos ip 7Ti (fc) + sin i\) 7T2 [k) , 
bk = — sin-i/; TTi(k) + cos-0 ^(k), 



(12) 



where 
cos 2-0 



CJfe - W T 



\/H - wt) 2 + 40 



(a;*. = ck/^/EE)- (13) 



With the new operators tt\ and 7T2, the Hamiltonian (11) 
is diagonalized: 



+ ^4 2, (^+(^ 2 (fc), 



(14) 



where the frequencies of the polaritonic branches are 
given by 



(*) = n 



uj k +WT±\/ (wife - wt) 2 + 40f 



(15) 



The Rabi splitting at resonance, i.e. at ujk = is 
20p. Eq. (15) allows to express the phenomenolog- 
ical parameter Op through the observables. Namely, 
Op = ^ujtujut/2 , where ujut "C is the transverse- 
longitudinal splitting. 

We note that the above description is valid only for 
wavenumbers k in the vicinity of "crossing" of the excita- 
tion branches, where u>k — rHT- It does not capture, how- 
ever, the correct behaviorlia of the polaritonic branches 
for k — > 0. In this limit an additional term of the type 
&kb-k + c.c. should be taken into account in Eq. (11). 



III. Diagonalization of the full Ha- 
miltonian 



where b is the annihilation operator of the medium ex- 
citations (e.g. excitons or optical phonons), which we 



Now let us consider the full Hamiltonian (1) with both 
Bragg scattering and light-medium interaction included, 
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H = J2 [w(«)a+ («)o^(g) + w(-g)a+ (-g)a^(-g)] 

+ o B X] t e# + «t(-9)^(g)] 

9 

+ cot^ [S+(«)S^(g) + 6+(-«)M-g) 

9 

+ o P X] [&i(«)fi^(«) + fii(«)M<?) 

9 

+ &t(-5)M-9) + at(-9)M-9) 



(16) 



where we have again truncated the "Bragg" Hamiltonian 
in (1) by including only near-resonance terms. If a col- 



umn of operators c = {a_>(q), a<-(—q), &-»(<?), b^(—q)} 
is introduced, then the Hamiltonian (16) can be formally 
rewritten in a matrix form 7i = c + Hc, where 



H 



( uj( q ) o B o P o \ 

B e-^ uj(-q) P 

P lo t 

\ o Op o uj T J 



(17) 



The four H eigenvalues yield the dispersion relations of 
the four excitation branches, whereas the eigenvectors 
determine the unitary transformation diagonalizing TL. 
The characteristic equation for the eigenvalues lo of the 
matrix H reads 



(uj(q) — lo^J (uj(—q) — u?j — 



ui(q) + u>(—q) — 2cj 



LOT — <*> 



of, 



{lot - w) 2 + Op = 0. 



(18) 



If the light-matter coupling is absent, i.e. Op = 0, then the roots of (18) reduce to two pure medium excitations 
with unperturbed frequency lot propagating in the forward and backward directions along z, and two purely photonic 
excitations with disperion relation given by Eq. (9) that results from the Bragg scattering. If, on the other hand, the 
Bragg scattering is absent, i.e. Ob = 0, then the roots of Eq. (18) reduce to two pairs of polariton branches with 
dispersion relation given by Eq. (15). 

It appears that the unitary transformation in four- dimensional space diagonalizing the Hamiltonian Tt can be 
parameterized by two angles: 



(19) 
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cos# sin-0 


sin 9 cosip 
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sin 9 cosip 
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9 sin?/; 


sin 9 simp 


— cos 9 cos tp 
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cos# simp 


sin# 


cos tp 


— cos cos ip 


— sin 9 sin ip 




X 




(-q)J 




\ sin 9 sin-i/j 


— COS 


9 cos?/' 


— sin 9 cos tp 


cos 9 sin?^ 







where Bi, D.2, and B2 are new operators that annihilate mixed light-matter states. The angle 9 in Eq. (19) is 
precisely the "Bragg" rotation angle introduced in Eq. (7) [for simplicity we set <p = for the modulation phase in 
Eq. (19)]. The second angle, ip is defined by the following relation 



cos 2ip 



yj [lo (q) - Lo(-q)f + 40| - 2(lo t - lo B ) 
|2(w T - Wb) - \l Hq) - u(-q)} 2 + 40| J +160 



(20) 



Naturally, for Ob = the angle ip reduces to the polari- 
tonic rotation angle "0 in Eq. (13). In the presence of the 
Bragg scattering, however, this rotation angle also de- 
pends on the "Bragg" parameters wb and B . Therefore 
it is the angle ip that characterizes the interplay between 
the polaritonic and diffraction effects. 



IV. Braggoritonic excitations 

In order to analyze the solutions of Eq. (18), it is conve- 
nient to introduce the following dimensionless variables. 
We measure frequencies, A, from the Bragg frequency, 
lob, and express them in units of the Bragg gap 20b: 



lob 



20 F 



(21) 



In analogy with Eq. (21) we introduce the dimensionless 
frequency detuning, S, of lot from the Bragg frequency 
lob, where 



S = 



LOT — 

20 R ' 



(22) 



As seen from Eq. (9), the natural unit for the wavevector 
deviation, q, from the Bragg wavevector, cr/2, is ctOb/wb- 
Hence, we introduce the dimensionless parameter 



(23) 



4 



.. 0.5- 

p 

u 

D 



p 
o 



-0.5 



(a) 5 = 

if 


"'' //V — «=° 

/ Sf a=0.35 

^ a=0.57 




8 





-1.5 -1 -0.5 0.5 1 

Dimensionless wavevector, Q 



1.5 




-1 -0.5 0.5 1 

Dimensionless wavevector, Q 

Fig. 1 Dispersion of mixed photonic-medium excitations for 
various coupling strengths: a = 0, 0.35, 0.57. (a) 5 = 0; (b) 
S = 0.3. The shaded area represents the two forbidden sub- 
gaps at a = 0.35. 

With the new notations, the excitation spectrum deter- 
mined by Eq. (18) can be rewritten in a more concise 
form 



\Q\ 




a" 



A-S 



(24) 



where a = J7p/(20b) characterizes the relative strength 
of the Bragg and polaritonic couplings. Expression (24) is 
our main result. It clearly demonstrates that the Bragg 
and polaritonic dispersion relations compete with each 
other. 

Consider for simplicity the case of exact resonance, 
i.e. S = 0. It is seen from Eq. (24) that in the absence 
of light-matter coupling (a = 0), the first term in the 
brackets gives rise to the conventional PBG. It is also 
seen that with increasing a (or, Op) the decay length, 
Im Q^ 1 , increases, and for sufficiently small A we find 
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Fig. 2 The effective mass for various excitations (in units 



of "free" Bragg mass Mg = 
ted vs. coupling strength: 



M 



Mg 2 at a 



solid lines are for 8 



0) is plot- 
= where 



Mgi = M B2 , = M} 2 ; dashed lines are for 6 = 0.3. 

that Q becomes real. This manifests the emergence of 
the novel allowed photonic states, or braggoritons, inside 
the PBG (see Fig. 1). The braggoriton branches in the 
excitation dispersion relations are described by the op- 
erators D.1 and D.2- They occupy the frequency ranges 
A 



0,±i(Vl + 16a 2 



1 



For small a (a <C 1), the 

braggoriton frequency interval reduces to ^0,±2a 2 ^. We 

note that due to the finite a value the Bragg gap broad- 
ens. Namely, the band edges of the branches described 
by the operators 'j&i , ify are respectively given for S = 

by A = ±\ (Vl + 16a 2 + l) (compare to A = ±1/2 for 

a = 0). The dispersion relations A(Q) calculated using 
Eq. (24) are shown in Fig. 1(a) for different values of a 
in the case of exact resonance wb = wt, or 5 = 0. 

Moderate frequency detuning 6^0 does not qualita- 
tively change the above picture as seen in Fig. 1(b). The 
major effect of frequency detuning is that the braggori- 
tonic branches 3i, Q 2 acquire an asymmetry since they 
are "pinned" by wt- Fig. 1 also shows that the Bragg- 
like photonic branches B\, B2 are affected by coupling 
or detuning only weakly. 

To quantitatively describe the braggoriton dispersion 
relation, we consider two characteristics: (i) dimension- 
less effective mass, M near the band edges that is defined 
from Eq. (24) by the relation, 



(25) 



and (ii) the density of states, N(A). Expanding Eq. (24) 
in Q yields the following effective masses for the brag- 
goritons 
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Fig. 3 Density of states for mixed photonic-medium exci- 
tations. Left panel is for the symmetric case (5 = 0), for 
a = 0.35 (solid line) and a = 0.57 (dashed line). Right panel 
is for 8 — 0.3 and a — 0.35. The thin solid vertical lines 
indicate the band edges. 
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^2 



\/(l-2<5) 2 + 16a 2 

1 + 25 
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(26) 



(27) 



These masses are plotted in Fig. 2 versus the coupling 
strength a. For a — > 0, we have M^, M^ 2 — > oo, re- 
flecting the fact that at a — the braggoritons reduce 
to dispersionless medium excitations that are not cou- 
pled to light. With the light-matter interaction switched 
on, the braggoriton effective mass rapidly decreases (the 
width of the in-gap branches increases). 

The one-dimensional density of states N(A) is given 
from Eq. (24) by: 



N(A) 



dA 



A - 




A-S 



(A - sy 



(28) 



and shown in Fig. 3 for different values of a and de- 
tuning, 6. The density of states of the braggoritonic 
branches inside the gap exhibits conventional ID square- 
root singularities at the band edges A = 5 and A = 



4a 2 



As mentioned above, the upper and lower Bragg-like 
photonic branches, B%, B2, are only slightly affected by 
the coupling and/or detuning. In particular their effec- 
tive masses, 



1 + 25 



16a 2 



(29) 
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(30) 



16a 2 



change only by a factor of 2 as a varies from zero to 
infinity (see Fig. 2) . 



V. Intragap localized states 

We now turn our attention to the localized photonic 
states caused by a phase-slip like defect. Note that in 
the absence of the polarizablc medium, a structure with 
one-dimensional modulation (3) of the dielectric function 
can be viewed as a distributed feedback resonatortl first 
considered by Kogelnik and |Shankc3 in 1972. Later it 
was realized that a phase slipE-l in the modulation 



e(z) = eo + Secos(az + <f)(z)), 



where 



(z) = 



61, z < 

62, z > 



(31) 



(32) 



results in a localized state inside the PBG. Within the 
second quantization formalism of Sec. II, the emergence 
of such a state can be established as follows. Consider the 
eigenstate annihilation operators of the Hamiltonian (4) 



01 

02 



cos 9 
sin 6 e~ 



sin 9 
cos 9 



(33) 



It follows from (33) that the absolute value, A, of the 
amplitude ratio of the left and right propagating waves 
constituting the eigenstates 0\, 02 is either A = tan#, or 
A = tan -1 9. These expressions are actually equivalent 
with appropriate choice of the sign of square root: 



A±(A) =2 A± 



(34) 



where we used the definition (7) of the rotation angle 9. 

In the presence of a phase slip (32) the continuity con- 
dition at z = reads 



A_(A)e-* = X*_(A)e 



(35) 



As is well knownEil Eq. (35) has a unique in-gap solution, 
A', for an arbitrary phase discontinuity <f>\ — 4>2, 



A' = cosx, 



where 



2 

h - 02 



-7T < 01 



<0, 



< 4>l — 4>2 < 7T. 



(36) 



(37) 
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Fig. 4 The frequencies of localized intragap states vs. the 
phase slip magnitude: Di, D2 are defect levels inside the 
Bragg gap (a = 0, dashed line); Di, D\ are defect levels in- 
side the lower subgap, and D2, Di are defect levels inside 
the upper subgap (5 — 0, a = 0.35, solid line). Thin solid 
and dashed lines represent the band edges of the two forbid- 
den gaps and of conventional Bragg gap (in the absence of 
coupling) , respectively. 

Generalization of the above consideration to include the 
polarizable medium is straightforward. It reduces to the 
following modification of the parameter A in Eq. (35): 



\±(A,a,6) 



2 I A — — - ± \ ( A ' 



A-5 



A-8 4 ' 



(38) 



Then condition (35) yields the gap state solution, A' 



& '=2 



5 + icosx± Ms - icosx 



4cv 2 



(39) 



where x is defined by Eq. (37). For a = we return 
to the in-gap state (36). Remarkably, we note that for 
nonzero coupling parameter a, when the Bragg gap is 
divided into two sub-gaps as in Fig. 1, the two values of 
A' determined by Eq. (39) are located in each of the cor- 
responding sub-gaps. In Fig. 4 the frequencies, A' of the 
localized in-gap states are shown versus the magnitude 
<pi — <f) 2 of the phase-slip for zero detuning <5 = 0. 

VI. Discussion 

Here we discuss the connection of the present work to the 
earlier related studies. 

The band structure of a photonic crystal with fre- 
quency dependent dielectric function was recently stud- 
ied numerically in Rcf. Eq, using the plane-wave method. 



In this work the photonic crystal was modeled as a two- 
dimensional array of GaAs rods. Frequency dispersion 
was introduced through the transverse- longitudinal split- 
ting of the optical phonons. The authorsc^l (see also Ref. 
pTj ) observed that numerous branches of the band struc- 
ture calculated for wb = become almost dispersionless 
at frequencies close to the frequency lut- 

In the context of the present work, this weakening 
of dispersion can be understood from Eqs. (26), (27) 
that describes the effective masses of the braggoritonic 
branches Ji, D.2- These masses rapidly increase as the 
coupling parameter a decreases. 

Signatures of braggoritonic excitations studied in de- 
tail in the present paper can be also found in the nu- 
merical calculations of Ref. [2^. In that work, transmis- 
sion spectra of a photonic crystal identical to that of 
Ref. |26J were calculated within the transfer-matrix for- 
malism. Two minima in the transmission spectra were 
found instead of the usual single minimum that is caused 
by the Bragg diffraction in non-dispersive photonic crys- 
tal. In light of the theory developed in the present work, 
these two minima can be identified with the two forbid- 
den subgaps in the excitation spectrum (Figs. 1 and 3). 
As we have demonstrated (Fig. 3), in the presence of 
light-matter coupling there are two spectral regions with 
zero density of states. Correspondingly, the transmission 
coefficient within these frequency regions must be low if 
the sample is sufficiently thick. 

A very different realization of periodic polarizable 
structures was the subject r©f|-extensive theoretical stud- 
ies during the last decade .r 9 n 3 ™ The structures are mul- 
tiple quantum wells separated by wide-gap semiconduc- 
tor barriers. The width, d, of each barrier was assumed 
to be close to Ao/2, where Ao is the wavelength corre- 
sponding to the intra-well exciton resonance frequency, 
wt- The condition d sa Ao/2 implies that the Bragg 
frequency, wb is close to wt- Since the quantum wells 
with strong frequency dispersion at w ~ lut had thick- 
ness much-smaller than d, then a real Bragg gap in the 
structuresEj'Ell was lacking. However, under the condi- 
tion lut = lu-q the dispersion law of light propagating 
along the principal axis was shown to have a gap within 
a frequency range \oj — ojt\ = (2r w T /7r) 1 ^ 2 = f2( eff ). 
Here r denotes the radiative rate for an exciton in a 
single well. In other words, il( cff ) plays, .the role of the 
"effective" Bragg gap in the structuresB3't2l Remarkably, 
a physical picture completely analogous to the multiple- 
quantum-well structures emerged from consideration of 
an optical lattice formed by laser-cooled atoms.Eil Corre- 
spondingly, the light dispersion relation derived in Ref. 
Bll has the same form as in 



n conclusion 
a detuning 



Refs. |29 

we note that as was recently pointed outl 
(lub — ^t) ~ f2( eff ) gives rise to a band of propagating 
states within the "effective" Bragg gap of the multiple- 
quantum- well structures. 
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